Decomposition of the Superwind in M82S 
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ABSTRACT 

We present new optical images {B, V, and Ha) of the archetypical 
starburst /superwind galaxy M82 obtained with the 8.2 m Subaru Telescope 
to reveal new detailed structures of the superwind-driven nebula and the 
high-latitude dark lanes. The emission-line nebula is decomposed into (1) a 
ridge-dominated component comprising numerous filament/loop sub-structures 
whose overall morphology appears as a pair of narrow cylinders, and (2) a diffuse 



^ Based on data collected at Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 
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component extended over much wider opening angle from the nucleus. We 
suggest that these two components have different origins. The ridge-dominated 
component appears as a pair of cylinders rather than a pair of cones. Since this 
morphological property is similar to that of hot plasma probed by soft X-ray, 
this component seems to surround the hot plasma. On the other hand, the 
diffuse component may arise from dust grains which scatter stellar light from 
the galaxy. Since inner region of this component is seen over the prominent 
"X" -shaped dark lanes streaming out from the nuclear region and they can be 
reproduced as a conical distribution of dust grains, there seems to be a dusty 
cold outflow as well as the hot one probed by soft X-ray and shock-excited 
optical emission lines. If this is the case, the presence of such high-latitude dust 
grains implies that neutral gaseous matter is also blown out during the course 
of the superwind activity. 

Subject headings: Galaxies: individual (M82) — Galaxies: intergalactic medium 
— Shock waves 



1. INTRODUCTION 

Bursts of massive star formation (starbursts) are considered to play important roles 
in the evolution of galaxies from high redshift to the present day. Approximately 10% of 
nearby galaxies show evidence for such starburst activity (Weedman et al. 1981; Balzano 
1983; Ho et al. 1997). In particular, some of them also show evidence for galactic-scale 
bipolar outflows (superwinds) driven by collective effects of a large number of supernova 
explosions occurred in the central region of galaxies (e.g., Chevalier & Clegg 1985; Tomisaka 
& Ikeuchi 1988; Heckman et al. 1990; Suchkov et al. 1994; Tenorio-Tagle & Munoz-Tunon 
1998; Strickland & Stevens 2000). The superwind phenomenon is also very important 
as well as the starburst activity itself because it provides metal-enriched gas into the 
intergalactic medium. 

Since the identification of optical emission-line filaments in M82 (Lynds & Sandage 
1963; Burbidge et al. 1964), much attention has been paid to this galaxy and it is now 
widely accepted that M82 is the archetypal starburst galaxy with evidence for the bipolar 
superwinds emanating from the central region to the outer halo area (e.g., Heckman et al. 
1987; Nakai et al. 1987; Bland & TuUy 1988; ShopbeU & Bland-Hawthorn 1998; Lehnert 
et al. 1999). Because of its proximity (a distance to M82, D — 3.63 Mpc; Preedman et 
al. 1994), M82 allows us to investigate detailed characteristics of the superwind activity. 
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Indeed, previous studies of the M82 superwind have shown that the gas in this superwind 
takes various phases: I) hot (T ~ lO''' K) soft X-ray emitting gas (e.g., Bregman, Schulman, 
& Tomisaka 1995; Strickland, Ponman, & Stevens 1997; Lehnert et al. 1999), II) warm 
(r ~ 10^ K) ionized gas probed by optical emission lines (e.g., Heckman et al. 1987; Bland 
&Tully 1988; McKcith et al. 1995; Shopbell & Bland-Hawthorn 1998), III) cool (T ~ 100 
K) neutral gas probed by dust grains, including i) dusty molecular-gas halo (e.g., Scarrott, 
Eaton, & Axon 1991; Sofue et al. 1992) and ii) complex dark lanes seen in the upper disk 
regions (Ichikawa et al. 1994; 1995), and IV) cold (T ~ 10 — 100 K) molecular and atomic 
gas (e.g., Nakai et al. 1987; Yun, Ho, & Lo 1993). However, it is not yet clear how such 
various gaseous components are related to each other although the observational properties 
of each component have been discussed extensively. 

In order to improve our understanding of this superwind activity, we present new 
sensitive high-resolution images of M82 taken with the 8.2 m Subaru telescope (Kaifu et 
al. 1998). These images are analyzed in various ways to reveal new detailed structures of 
the galactic disk (stars), dark lanes (dust grains), ridge-dominated and diffuse components 
(shock-heated warm gas and dust grains scattering the nuclear light, respectively) of the 
nebula. Combining with the previous information of X-ray image for hot ionized gas, we 
investigate the multi-structured superwind phenomenon in M82. 

2. OBSERVATIONS AND DATA ANALYSIS 

The optical imaging observations of M82 were made with the FOCAS (Faint Object 
Camera and Spectrograph: Kashikawa et al. 2000, Yoshida et al. 2000) on the Subaru 
Telescope during a commissioning run on 2000 February 2 (UT). The seeing condition was 
0.7" - 0.8" for all images. In order to investigate the emission-line nebula associated with 
the superwind, we obtained direct images with two narrowband filters, N658 and N642. 
The N658 filter centered at 6588A with a band width of 73A was used to probe both Ha 
A6563 and [N ll] AA6548, 6583 emission lines (hereafter. Ha + [N ll]), while the N642 filter 
centered at 6428A with a band width of 127A was used to measure the adjacent continuum 
emission. We obtained five N658 images and three N642 ones. Each exposure time was set 
to be 120 seconds and thus the total integration was 600 seconds and 360 seconds for the 
N658 and N642 images, respectively. We also obtained optical B- and y-band images. The 
total integration time was 120 seconds (4 x 30 seconds) and 75 seconds (3 x 25 seconds) 
for B and V, respectively. No photometric standard star was observed and thus we do not 
calibrate the fiux scale of any images. All the data reduction was made using the dedicated 
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software for FOCAS implemented with IDL (Yoshida et ah 2000) and IRAFi. 

These images were further processed in the following ways. A pure Ha + [N ll] 
emission-line image (hereafter Ha + [N ll] image) was made by subtracting the scaled 
N642 image from the N658 one i . A pure continuum image around 6500A (hereafter 
6500A continuum image) was made by summing N642 and N658 images and subtracting 
the Ha + [N ll] image from the summed one. The Ha + [N ll] image is shown superimposed 
on the true-color images of the M82 galaxy disk (made from B, V, and 6500A continuum 
images) in Figure 1. The emission-line nebula appears to consist of the following two major 
components; (1) the ridge component which corresponds to the bipolar filament-dominated 
component, and (2) the diffuse component which is more widely extended around the ridge 
component. 

In the ridge component (Figure 2), many shell- and loop-like filaments are recognized. 
Their typical sizes are 0.8" (unresolved) - 2" wide and 5" - 20" long. In order to show 
such filaments more clearly, we developed a ridge detection program which searches for 
any regions with large negative value of the second derivative of the image (i.e., tracing 
local peaks in the emission- line image). The sky area without filaments was extracted from 
the Ha + [N ll] image by masking out the ridges detected, and the image of the diffuse 
component was made by filling-in the masked regions by interpolating the surrounding 
pixel values. Subtracting the diffuse component image from the Ha -|- [N ll] image, we 
obtain an image of the ridge component. Figure 3 shows images of Ha + [N ll] nebula 
shown separately in ridge and diffuse components. Note that since our ridge-detection 
program can only detect the extreme limbs or tangent points of the structures, structures 
resolved slightly around the peak of the ridges are included in the diffuse component, and 
are notable especially near the root of the nebula. Therefore the amount of the diffuse 
component is overestimated. The broadband color image was made by dividing the i?-band 
image with the 6500A continuum image, and is shown in Figure 4 along with the ridge and 
diffuse components of the nebula. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

■^The scaling of the N642 image was done by comparing fluxes of nearby stars around M82. Although we 
did not take account of the possible underestimate of the scaled N642 flux due to Ha absorption in N658 
image, this effect is estimated to be negligibly small since the disk stellar emission could be subtracted clearly 
over the entire disk. 
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3. RESULTS and DISCUSSION 

3.1. The Ridge Component 

The ridge component comprises many filaments along the superwind. Each filament 
shows either shell- or loop-like structures directing radially from the nuclear region. Similar 
filament-dominated emission-line nebulae are also found in some well-known Galactic 

objects with an outflow [such as Ml-67: a Wolf-Rayet star with stellar wind (Grosdidier et 
al. 2001) and M57: a planetary nebula with a ring (Komiyama et al. 2000)]. It is widely 
accepted that such a filament can be a local shock front between the inner expanding matter 
and the outer interstellar medium. The observed filaments in M82 are also considered to be 
driven by shocks. This interpretation is supported by the optical emission-line diagnostics 
(Heckman et al. 1987) although their data probe both ridge and diffuse components. 

As mentioned before, the overall structure of the ridge component looks like a pair of 
elongated cylinders rather than a pair of cones, and hence hereafter we simply denote it as 
the cylinder. The diameter of the cylinder at the galaxy disk plane (35" - 40") is similar 
to, or slightly larger than, that of the nuclear starburst region traced both by near-infrared 
emission (30" - 35"; e.g., Lester et al. 1990) and by an ensemble of non-thermal point-like 
sources (e.g., Kronberg, Biermann, & Schwab 1985). It seems worthwhile noting that the 
cylinder shows the sharp boundary, in particular at NE, SE, and SSE of the nebula, causing 
the limb-enhanced morphology. All these features can be well reproduced by numerical 
simulations of the superwind in which radiative shock occurs around the expanding hot gas 
(e.g., Suchkov et al 1994; Strickland & Stevens 2000) (see section 3.3). 

The hot gaseous component probed by soft X-ray emission is known to be coUimated 
rather tightly along the minor axis of the galaxy than the optical emission-line nebula 
(Bregman et al. 1995; Strickland et al. 1997; Lehnert et al. 1996). We found that the ridge 
component surrounds the hot gaseous nebula more closely than the diffuse component. 
However, the tangential profiles along the galaxy disk at high latitude look different between 
the ridge component (showing limb-enhanced profile) and the soft X-ray nebula (showing 
the ridge along the centerline). Note, however, that such difference can be explained if 
we assume that the optical filaments are distributed on the thin surface around the inner 
hot gas component. All these lines of evidence suggest that the radiative shock at the 
thin interface between the expanding hot X-ray emitting gas and the ambient cold gaseous 
matter would be responsible for the ridge component. 
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3.2. The Diffuse Component 

The diffuse component was first found by Bland & Tully (1988) and was confirmed by 
Shopbell & Bland-Hawthorn (1998). It shows rather oval shape elongated along the minor 
axis of the galaxy as a whole, but is divided into NW and SE parts around the disk plane. 
It is extended more widely in the tangential direction (parallel to the galaxy disk) and does 
not show a limb-enhanced morphology, being much different from the ridge component @ . 
These observational properties suggest that this component does not have a form of thin 
surface but occupies a large part of the volume probed by optical emission lines. 

Here a question arises as what the origin of the diffuse component is. Previous imaging 
polarimetries of the M82 nebula have revealed strongly-polarized emission whose E vectors 
show circular symmetric pattern around the nucleus. Since this property is a typical 
signature of the reflection nebulae, the dust scattering of the galaxy light would contribute 
to the nebula (e.g., Bingham et al. 1976; Schmidt, Angel, & Cromwell 1976). It is thus 
suggested that the diffuse component arises from the scattering of the galaxy light at dusty 
halo. 

In addition to the dusty halo, M82 is well known to have rich dark lanes across the 
galaxy disk. They are peculiar and complex in shape, and are remarkably different from 
those found in normal spiral galaxies. The most prominent and peculiar dark lane is seen 
on SE side of the disk (the front side of the disk: e.g., Shopbell & Bland-Hawthron 1998; 
McKeith et al. 1995). It runs through the nuclear region out to 2.5 arcmin (2.6 kpc) from 
the nucleus (Figure 5). Interestingly it curves toward both SE and SW directions at both 
ends of the dark lane, and can be traced up to ~ 1' (~ 1 kpc in projection) above/below 
the galaxy disk at their tips. Similar but less-prominent dark lanes are also found at another 
side of the disk emanating to NE and NW from the nuclear region. These two sets of the 
dark lanes make an "X" -shaped morphology around the nucleus. They are considered to be 
distributed not within the galaxy disk but above/below the disk (e.g., Ichikawa et al. 1994; 
1995). We find that the diffuse component can be seen even on the dark lanes, suggesting 
that dust grains responsible for the dark lanes are also responsible for the scattering of the 
diffuse component. 



''As noted in the earlier section, separating the ridge and diffuse components by the ridge-detection 
program could not be perfectly done because some part of the flux separated as "diffuse" may actually be 
structures which arc resolved slightly around the peak of the ridges. However, since the ridge component 
sppares to be a collimated cylindrical structure outer parts of the nebula around the ridge component are 
likely to be truly "diffuse" . Therefore the comparison of the overall shapes between the ridge and diffuse 
components can be made with less ambiguity. 
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The next question is addressed to the origin of the dark lanes. The dark lanes have 
the following characteristics. First, the X-shaped dark lanes seem to be streaming out 
from the nucleus. Second, the NW boundary (and also probably the NE boundary) of the 
SE dark lane is smoothly curved comparing with other complicated dark lanes. Finally 
one-dimensional profile of the SW dark lane, cut perpendicular to the disk across the 
dark lane (Figure 5), shows highly asymmetric structure of the extinction, i.e., the dust 
extinction increases monotonically toward the NW side, and then decreases suddenly at the 
NW edge the SE dark lane. Such characteristics can be understood if we assume that the 
dust is distributed in a wide cone-like geometry and we see a higher extinction near the 
cone surface due to projection effect. This supports a scenario of the dusty superwind in 
which dust grains arc provided from the nuclear region with the outflow associated with 
the superwind probed in optical emission-lines and soft X-ray emission. If this is the case, 
high-latitude dusty halo may be created in a similar way as for the dark lanes (i.e., the 
dusty superwind) since there arc no distinct boundaries in the diffuse component dividing 
the dust scattering on the dark lane near galaxy plane and high-latitude halo area up to 
^ 2 kpc from the disk. Yun et al. (1993) suggested that a part of the high-latitude HI 
gas of M82 may be affected dynamically by the superwind, supporting for the presence of 
a high-latitude cold gas associated with the superwind outflow. Although we do not rule 
out an alternative possibility that the past tidal interaction with the neighbor galaxy M81 
has brought some dusty matter from the galaxy disk up to the halo (Yun, Ho, & Lo 1993; 
1994), the observational properties discussed in this paper suggest that the dusty superwind 
scenario is more reasonable. 

Such a dusty cold-gas outflow can be made as a result of the pushing out and/or 
dragging out process of the dusty gas around the nucleus by the expanding hot high-pressure 
gas within the superwind (e.g., Suchkov et al. 1994; Strickland & Stevens 2000). Indeed 
dusty outflows associated with supcrwinds have been suggested by optical observations 
(Nakai 1989; Heckman et al. 2000; Veilleux & Rupke 2002; Rupke, Veilleux, & Sanders 
2002) and submiUimeter observations (Alton, Davies, & Bianchi 1999). 

Finally we give a comment on the previous claim that even the inner fllamentary 
nebula shows strong polarization whose direction is nearly aligned to that of the outer 
diffuse nebula (Scarrott et al. 1991), which could be inconsistent with our conclusion 
since no polarized flux is expected from the shock-excited emission. We suspect that 
the imaging polarimetry with insufficient spatial resolution could be the cause of this 
apparent contradiction. Separate polarimetric measurement of the filaments in the ridge 
and diffuse components are required to disentangle the polarimetric properties of the narrow 
filamentary structures (typically ^ 2" in width) and the overlapping diffuse component. 
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3.3. Implications for Superwind Model 

Important characteristics of the M82 superwind nebula are that the shock-driven 
ridge component contains numerous shells and filaments and its overall shape is rather 
coUimated out to outer region of the nebula. In recent hydrodynamical simulations of 
the superwind conducted by Strickland & Stevens (2000), two types of the interstellar 
medium (ISM) density distribution, thin or thick disk, are adopted to see its effect on 
the evolution of the wind. It turns out that the model with the thick ISM can reproduce 
the characteristics of M82 nebula: The thick and dense wall of ISM created around the 
starbursting region works to coUimate the wind, and Rayleigh- Taylor instability at the 
wall creates small-scale structures which are eventually dragged out with the wind to form 
shells/filaments (Strickland & Stevens 2000). Although these hydrodynamical simulations 
with thick ISM look successful, we point out the following points which require some 
cautions. As Strickland & Stevens (2000) noted, the thick ISM cannot reproduce the 
observed rotation curve of M82. Therefore one may need more realistic ISM distribution 
models or improved treatment of gas dynamics at the circumnuclear region to reproduce 
both a coUimated wind and the rotation curve consistently. It seems also noteworthy 
that the morphology of the small-scale structures of the nebula depend rather strongly 
on the numerical resolution in the simulation, and higher resolution tends to result in 
finer shell/filament structures (Strickland & Stevens 2000). Therefore one may need even 
finer-scale numerical simulation to enable us to compare sizes and spatial distributions 
of shells/filaments between observations and the models. Note that, although observed 
small-scale structures with a typical size of 90 — 350 pc (5" - 20") are already resolved by 
their simulations with a cell size of 4.9 — 14.6 pc, the simulation resolution would not be 
fine enough to resolve the circumnuclear region where smaller-scale instability begins to 
work. Finally, we show a cartoon in which all the superwind components discussed here are 
summarized (Figure 6). 

We would hke to thank all staffs of the Subaru telescope. 
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Fig. 1. — True color image of M82 made from B, V, and 6500A continuum images. The total 
Ha + [N ll] image (ridge plus diffuse components) is shown superimposed. Some horizontal 
lines and circles are the artifacts in the image processing. See a main text for the detail of 
the image processing. 

Fig. 2. — Some interesting small-scale structures within southeastern superwind nebula. 
Central panel shows Ha + [N ll] image of the ridge component. The ridge component image 
is rotated 65 degree clockwise so that southeastern nebula is seen in horizontal direction and 
the stellar disk is seen at vertical direction at the right edge of the panel. Selected regions 
are marked with blue boxes, and each of them is labeled from A to F. Six sub-images around 
the central one (labeled as A - F) show enlarged images of each corresponding sub-region 
marked in the central panel. Total (ridge plus diffuse) component is shown in these panels to 
present full detail of the structures. All images are shown in a same scale, and their box size 
is 30 X 30 arcsec^ for all images except for the panel F whose size is 30 x 60 arcsec^. All images 
(central one and six sub-images) are shown in square-root intensity scale (as shown in lower 
right color bar), and maximum and minimum intensity cut levels are adjusted to present the 
fine structures clearly depending on the positions and the intensity of each structure. 

Fig. 3. — Ha -|- [N ll] images decomposed into ridge and diffuse components. Upper three 
panels show entire nebula (4.0 x 5.0 arcmin^) of ridge component (left), diffuse component 
(middle), and total (ridge plus diffuse) component (right). These three images in the upper 
row are shown in logarithmic scale. Lower three panels show inner 1.5 x 2.0 arcmin^ region 
of the upper image. Ridge, diffuse, and total components are shown separately as in the 
same way as for the upper panels. These three images in the lower panels are shown in 
square-root scale. Tick marks are shown every 10" for all images. 

Fig. 4. — Comparison of various components of at central 3.5 x 3.5 arcmin^ regions of M82. 
(upper left) the 6500A continuum image representing the stellar distribution along the 
galaxy disk shown in square root scale, (upper right) the S-band image divided by the 
6500A continuum image representing the distribution of dark lanes shown in linear scale, 
(lower left) the ridge component of the emission-line nebula representing the shock-driven 
nebula shown in logarithmic scale, (lower right) the diffuse component of the emission- 
line nebula representing the dust scattering component of the nebula shown in logarithmic 
scale. An outline of the X-shaped dark lanes (upper right) is superimposed on the diffuse 
component. Tick marks are shown every 10" for these four images. 

Fig. 5. — One dimensional profile of the dark lane. The cut region runs through I'SWW of 
the nucleus along the disk, and extendes over 1.5' perpendicular to the disk. The profile 
is averaged over 10" along the disk. The x-axis is shown in a relative distance from the 
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mid-plane of the disk in arcseconds, and y-axis is shown in a B-{red continuum) color in 
magnitude with a tick corresponding to 0.1 mag. See text for the red continuum. 

Fig. 6. — Schematic view of M82 seen from the side. 
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